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M.  ABSTRACT  fCanMaua  0 rararaa  aMa  H nocoaaary  and  t dam  tty  by  block  numbat ) 

A new  class  of  resonators  has  recently  been  proposed  which  has  the  pro- 
perty of  having  different  magnifications  M and  M in  orthogonal  directions 

x y 

perpendicular  to  the  propagation  direction;  the  most  interesting  cases  are 
those  with  collimated  output  analagous  to  conventional  confocal  unstable  reso- 
nators. The  designs  utilize  spherical  mirrors  which  are  tilted  substantially 


(ABSTRACT  continued) 


DD 


FORM 

I JAM  71 


1473 


EDITION  OF  I MOV  6S  IS  OBSOLETE 


CLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Bhan  Data  Bntarad) 


1 


J 


SECURITY  CLASSIFICATION  OF  THU  PAOEfWhan  Data  Bnttnd) 


(ABSTRACT  concluded) 

from  the  mirror  normal.  The  conventional  orientation  requires  the  mirror  normal 
to  be  along  (or  very  nearly  parallel  to)  the  propagation  direction.  Use  of 
spherical  mirrors  by  this  new  class  of  resonators  eliminates  the  difficulties  of 
fabricating  mirrors  with  more  complex  figures,  which  would  otherwise  be  required 
to  achieve  the  desired  asymmetric  magnification.  The  critical  problem  of  keeping 
beam-quality  degradation  (aberrations)  due  to  mirror  tilt  within  acceptable  limits 
has  been  examined  experimentally  and  analytically.  The  results  confirm  the 
possibility  of  achieving  acceptable  small  aberrations  for  large  tilt  angles 


(convex  mirror  tilt  as  large  as  70< 
operation  at  infrared  wavelengths. t 


, corresponding  to  140°  beam  deflection)  to 

-dec. 


A confirmation  of  the  tilted  spherical  mirror  resonator  concept  was  val- 
idated theoretically  and  experimentally  by  choosing  a convex/ concave  mirror  set 
with  small  radii  of  curvature  (small  f number) . The  deliberately  exaggerated 
aberrations  were  compared  by  experimental  interferometry  and  computational  plots 
of  optical  path  differences.  The  results  of  experiment  and  theory  agreed  quan- 
titatively. — . 

The  following  conclusions  were  drawn  from  this  work: 

1)  The  sense  of  relative  tilt  of  the  spherical  mirrors  affects  the  aberra- 
tions. 

3 

2)  The  dominant  terms  contributing  to  aberrations  are  of  the  form  y and 

2 3 

x y with  y predominating  for  high  aspect  ratio  resonators. 

3)  The  concept  has  many  advantages  over  toroidal  resonators  used  with 
linear  chemical  lasers.  Chemical  laser  development  programs  experiment  with 
varying  values  of  x^,  thereby  requiring  several  sets  of  toroidal  mirrors  for 

beam  extraction.  The  present  concept  has  the  ability  to  vary  M continuously 
to  account  for  planned  variations  in  x of  interest  in  chemical  laser  develop- 
ment experiments,  thereby  requiring  no  more  than  one  resonator. 
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This  new  class  of  resonators,  which  possesses  certain  unique  and 

1 2 

desirable  features,  was  proposed  and  discussed  earlier.  ’ The  basic 
design  equations  were  derived  and  experimentally  confirmed,  and  a report 
was  given  of  a preliminary  interferometric  study  of  optical  quality 
expected  from  such  resonators.  This  class  of  resonators  is  useful 
provided  only  that  the  output  optical  quality  can  be  maintained  at  an 
acceptably  high  level.  A more  detailed,  combined  theoretical-computa- 
tional study  of  expected  optical  properties  has  subsequently  been  per- 
formed. It  forms  the  subject  of  this  report. 

This  study  considers  the  effect  of  one  complete  passage  from  convex 
to  concave  mirror  on  an  initially  collimated  beam.  This  is  essentially 
the  conventionally  aligned  travelling-wave  (or  ring)  resonator  configu- 
ration. Standing-wave  resonators,  with  or  without  beam  inverters,  are 
also  of  considerable  interest,  but  they  are  not  considered  here.  Both 
cases  of  relative  sign  of  spherical-mirror  tilt  angles  are  considered; 
these  are  referred  to  as  U and  Z configurations.  The  Z configuration 
is  preferred  over  the  U configuration  in  optical  quality. 

A two-dimensional  treatment  (i.e.,  one  which  neglects  one  of  the 
two  transverse  dimensions)  is  useful  as  a conceptual  vehicle  and  as  a 
means  of  obtaining  useful,  though  only  approximate,  estimates  of  output 
beam  quality.  This  is  discussed  first.  A three-dimensional  ray-tracing 
approach  is  then  discussed  and  its  predictions  are  compared  with  experi- 
ments. The  most  important  part  of  this  study  is  the  experimental  study, 
which  is  an  extension  and  refinement  of  the  preliminary  interferometric 
study  reported  earlier. 


II.  THEORETICAL  CALCULATIONS 

A.  Two-Dimensional  Approximation 

The  two-dimensional  approximation  has  the  advantages  of 
greater  ease  of  detailed  calculations,  easier  visualization,  and  feasi- 
bility of  obtaining  approximate  closed-form  predictions  of  optical 
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quality.  On  the  other  hand,  it  is  certainly  not  highly  accurate  quanti- 
tatively, and  it  sheds  no  light  on  the  appearance  of  interferograms  and 
their  behavior  with  respect  to  variation  of  experimental  parameters. 


An  approach  based  on  local  wave-front  curvature  was  pursued  first; 
this  led  to  a closed-form  approximate  prediction  of  the  magnitude  of  the 

3 

principal  contribution  (the  y term)  to  the  aberration  resulting  from 
tilting  the  spherical  mirrors.  The  assumptions  and  approximations 
underlying  this  approach  are  perhaps  obscure,  so  calculations  were  also 
carried  out  for  one  particular  set  of  parameters  using  a ray-slope  and 
an  optical-path-difference  (OPD)  approach.  The  latter  two  calculations 
agreed  with  results  of  the  former  (local  wave-front-curvature)  approach, 
thereby  lending  it  credibility. 

The  two  basic  types  of  configurations,  U-type  and  Z-type,  are 
illustrated  in  Figures  1(a)  and  1(b).  The  basic  design  equations 
derived  earlier  make  no  distinction  between  these  two  types  of  design. 

The  various  treatments  of  aberrations  as  given  here,  however,  clearly 
distinguish  between  the  two  types. 

1.  Wave-Front  Curvature  Approach.  Figures  1(a)  and  1(b) 
present  the  two  basic  types  of  configurations.  The  combination  of 
mirror  curvatures,  separation,  and  tilt  angles  for  the  central  ray  will 
be  chosen  such  that  the  curvatures  in  the  y-z  and  x-z  plane  introduced 
at  the  convex  mirror,  after  modification  by  passage  to  the  concave 
mirror,  (will  be  compensated  by  the  concave  mirror),  and  the  emergent 
central  portion  of  the  beam  will  again  be  collimated.  This  condition 
will  not  be  exactly  satisfied  for  noncentral  portions  of  the  beam,  as 
schematically  related  to  off-center  rays  sketched  in  Figures  1(a)  and  1(b). 
Behavior  at  noncentral  portions  of  the  beam  between  curved  mirrors 
differs  from  that  at  the  central  position  in  two  respects: 

a)  The  local  tilt  angles  of  the  mirrors  are  slightly  different 
from  the  central  values. 

b)  The  distance  between  points  of  reflection  by  the  two  mirrors 
is  slightly  different  for  the  off-center  case. 

The  increments  in  tilt  angles  at  the  convex  and  concave  mirrors, 
respectively,  in  going  to  an  off-center  case  are  labelled  60^  and  602 

initially.  These  two  values  are  (approximately,  but  to  a high  order) 
equal.  This  can  be  shown  as  follows: 
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From  Equation  (32)  , it  is  noted 


By  Equation  (28)  , it  follows  that 


sec  0, 


sec  0, 


Combining  these  two  equations  gives 


The  positive  sign  applies  to  the  U configuration  (from  its  shape);  the 
negative  sign  applies  to  the  Z configuration.  For  simplicity,  6 is  used 
to  denote  the  common  magnitude  of  60^  and  602* 

a.  Increment  in  Separation  Required  to  Maintain  Con- 
focality:  The  basic  confocality  condition  required  to  obtain  collimation 

of  the  output  beam  in  a localized  region  is 


L (6=0) 
req 


J [-  lRil 


cos  0^  + R2  cos  ®2 


The  incremental  value  [i.e.,  change  in  L required  to  maintain  confoca- 
lity as  the  ray  of  interest  moves  from  the  center  toward  the  upper 


Op.  cit . , Cason  et  al..  Optical  Resonators 


Ibid. 
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to 


portion  of  Figures  1 (a)  and  1 (b) ] can  be  obtained  by  expanding  about 
the  central  angle.  For  the  incremental  ray,  60^  is  considered  the 

negative  of  6 (i.e.,  6 is  a positive  quantity).  This  gives 


4Lreq  ’ 1 6 [*  lRlJ  el"  h - R2  sIn  h] 


b.  Actual  Increment  in  Separation  of  Reflection  Points: 
First,  only  the  difference  in  Z-coordinates  of  the  points  of  reflection 
of  a ray  from  each  of  the  mirrors  is  considered.  By  geometrical  con- 
struction, it  can  be  seen  that 


AZ  »*  - 6^-  | R1 1 sin  6^  + R2  sin  0^  ^ 


From  separate  numerical  examination,  it  is  found  that  the  increment 
in  path  length  due  to  the  change  in  y coordinates  is  small  compared  to 
the  preceding  change  in  z coordinate.  The  incremental  focal  length 
then  is  the  difference  between  these  quantities,  and  it  follows  that 


AF  = AL  - AZ 
req 


I 6 [-  !Ri! 


sin  0^  + R^  sin 


c.  Radius  of  Curvature  of  Emerging  Wavelet:  Next  the 
local  radius  of  curvature  of  the  wave  reflected  from  the  concave  mirror 
at  the  displaced  position  is  determined  for  a mirror  spacing  adjusted 
such  that  the  center  portion  of  the  emergent  wave  is  collimated.  From 
the  lens  equation,  it  follows  that 


<yaF) 


<V 


where  F^^  is  the  y-coordinate  effective  focal  length  of  Mirror  2 (con- 
cave mirror);  namely,  1/2  R?  cos  Q^. 

If  it  is  assumed  that  the  equation  for  the  curve  of  constant  phase 
(obtained  by  setting  x = 0 in  the  equation  for  the  surface  of  constant 
optical  phase)  is 


z - A y = 0 , 

the  following  is  readily  obtained: 


= 6 A y 
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By  definition,  the  left-hand  side  of  this  equation  is  the  reciprocal  of 
the  local  radius  of  curvature.  That  is 


h = 6 A y 
R0 

Substitution  in  the  preceding  equation  gives 


phase.  Introducing  the  earlier  expression  for 


where  AF  is  given  by  the  preceding  expression.  Here  y is  used  to  denote 
the  displacement  at  the  output  position;  hence,  y = 6c os  e Com- 

bining all  of  these  gives 

z = -3  V"  [-  \\\  sin  9 ± R sin  ejy3 

R2  cos  ®2  U J 

By  inspection,  the  value  of  2 for  a given  value  of  y is  always  smaller 
in  magnitude  for  the  Z configuration  than  for  the  U configuration. 

(Only  the  magnitude,  but  not  the  sign,  of  z as  predicted  by  this  equa- 
tion is  of  significance,  because  the  value  of  z has  opposite  signs  for 
rays  incident  on  opposite  sides  of  the  convex  mirror) . 

As  a numerical  example,  a case  with  R^  = -290  cm,  R^  = 675  cm, 

= 55°,  y = 9.176  cm  (corresponds  to  6 0.01)  is  considered.  The  asso- 
ciated value  of  is  38.76°.  Using  the  preceding  formula  for  the  Z 
configuration,  the  following  is  obtained: 

z = 1 x 185.035  x (9. 176)  3 = 0.000  980  cm 

(675  x 0.77977) 

For  the  U configuration,  the  following  is  obtained 

Z = (675  x~  0.77977)  * (~660-i44)  * (9. 176) 3 = -0.003  497  cm 


for  the  curve  of  constant 
Rq  gives 


AF  2 

- y 


3 n2  2 „ 
r2  cos  e2 
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The  expression  for  z as  a function  of  y can  alternatively  be  written  as 
an  explicit  function  of  the  magnifications  and  , rather  than  the 

angles  e and  q^,  although  this  has  not  been  done  here.  There  are 

equations  from  which  0 can  be  determined,  when  0.,  R.  , and  R„  are 
^ 2 1 i z 

specified . 


This  detailed  discussion  has  been  given  for  the  wave-front  curva- 
ture approach  in  the  y-z  plane  because  it  gives  an  explicit  expression 
in  terms  of  resonator  parameters.  Alternative  approaches  which  have 
been  used  will  be  more  briefly  discussed.  The  function  of  the  latter 
has  been  essentially  limited  to  numerical  and  analytical  confirmation 
of  the  preceding  approach  for  one  set  of  resonator  parameters. 

2.  Ray  Slope  Approach.  A two-dimensional  ray-trace  calcu- 
lation was  carried  through  numerically  for  the  reference  case  with 
R^  = -290  cm,  R^  = 675  cm,  0^  = 55°;  here  the  small  parameter,  again 

labelled  6 (the  distance  between  points  of  intersection  of  two  incident 
rays  with  convex  mirror  and  measured  in  units  of  radius  of  curvature  of 
convex  mirror)  was  treated.  Numerical  results  agreed  with  the  wave- 
curvature  results.  Treatment  was  done  semi-analytically  with  various 
quantities  expressed  as  expansions  in  powers  of  6.  From  the  latter, 

3 2 

the  y dependence  of  z (y  dependence  of  slope)  was  confirmed  as  found 
from  the  wave-curvature  approach. 


3.  OPD  Approach.  Direct  numerical  and  semi-analytical 
expansion-type  calculations  were  made  to  determine  the  OPD  for  the 
previously  mentioned  reference  case.  Results  again  confirmed  the  wave- 
curvature  results. 

4.  Beam-Quality  Predictions.  Deterioration  of  the  optical 
quality  of  the  output  beam  is  essentially  determined  by  the  root-mean- 
square  (rms)  value  of  variations  in  optical  path  length,  z,  measured  in 
units  of  wavelength,  relative  to  a suitable  reference  plane  for  a colli- 
mated output  beam.  The  rms  variations  are  considerably  smaller  than 
peak-to-peak  or  center-to-peak  variations,  and  the  choice  of  reference 
plane  is  important. 


z 

cp 


In  the  two-dimensional  approximation,  the  center-to-peak  variation, 
(approximately  half  the  peak-to-peak  variation)  can  be  expressed  as 


z 

cp 


03 


J03 


y 
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^Ibid. 
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is  a certain 


A *- 


where  y is  the  half-width  of  the  output  beam,  and  C„_ 
max  03 


expansion  coefficient  (use  of  the  double  subscript  is  for  consistency 
with  the  later  three-dimensional  calculations) . From  the  preceding 


equation,  the  approximate  prediction  for  is 


C03* 


r3  cos3  e2 


sxn 


+ R2  sin  ©2 


The  rms  value  for  an  untilted  reference  plane  is  obtained  by  simple 


integration 

2 


y 

f 

-y 


[z  (y)]' 


dy 


rms,nt 


2y 


03 

2y 


J € 

-y  y 


dy 


rms.nt 


1 r 2 6 
7 C03  y 


The  factor  of  7 reduction  is  naturally  of  considerable  practical  sign- 
ificance. It  follows  that 


rms , n t 


0.378  z 


cp 


The  preceding  equations  assume  a reference  plane  which  is  perpen- 
dicular to  the  central  ray  of  the  output  beam.  There  is  no  a priori 
reason  why  the  reference  plane  must  be  so  chosen.  A reference  plane 
which  optimizes  the  beam  quality  can  be  chosen  (i.e.,  to  allow  for  a 
tilted  reference  plane  (this  is  closely  related  to  allowance  for  a 
tilted  reference  beam  in  interferometry,  as  discussed  later),  as  given 
by  the  equation 


z , = C.  y 
ref  t 


The  integral  involved  is 


rms 


,t  * i i,y  [co3  y3  - v] 

' k C032  4 [y  l - B>j‘  dy 


dy 
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where 


B = 


J03 


The  following  is  readily  obtained: 


‘rms.L2  ' 7 C032  y6  £<B) 


where 


f (B)  = 1 - B„  + 7 B 


5 y2  3 yA 


To  minimize  f(B),  its  derivative  is  taken  with  respect  to  B and  set 
equal  to  zero.  The  following  is  obtained: 

B .1**  • 

opt  5 


The  minimum  value  of  f (B)  is  then 
fopt<B)  -2-5  ■ °-160 

This  represents  a further  substantial  reduction  in  the  OPD  which  is 
effective  in  deterioration  of  beam  quality.  For  optimal  reference- 
plane  tilt,  the  following  is  obtained: 


'rms, opt 


-A.  c 2 y6 
175  03  y 


The  minimized  rms  value  itself  can  then  be  written  as 


z 


rms , opt 


= 0.151  C 


03 


3 

y 


This  is  less  by  a factor  of  approximately  13  than  the  peak-to-peak  OPD. 
It  is  the  square  of  this  factor  which  is  important  in  determining  beam- 
quality  deterioration. 

Some  numerical  examples  are  given  in  Tables  1 through  3,  which  are 
representative  of  the  mirror  pair  studied  interferometrically : an 
example  chemical-laser  resonator  design,  and  a CO^  pulsed  electron-beam 


laser,  respectively.  From  these  estimates,  the  optical  quality  should  be 
quite  acceptable  at  infrared  wavelengths,  when  rms  and  reference-plane- 
tilt  effects  are  accounted  for,  except  for  the  rather  extreme  case  of 
the  mirrors  being  used  for  the  interferometric  testing  in  the  case  of 
Table  1.  This  pair  of  mirrors  is  considered  at  full  aperture  primarily 
because  they  provide  useful  interferometric  structure  at  visible  wave- 
lengths, thereby,  facilitating  quantitative  experimental  studies  of 
aberrations . 


TABLE  1.  TWO-DIMENSIONAL  APPROXIMATE  OPD  PREDICTIONS  (IN  MICRONS) 
FOR  A Z CONFIGURATION  TRAVELLING-WAVE  RESONATOR  WITH  K1  = -290  cm, 

R„  = 675  cm,  Y = y =7.62  cm,  FOR  VARIOUS  VALUES  OF  CONVEX- 
2 max 

MIRROR  TILT  ANGLE  Gj* 


01  (deg) 

z 

cp 

z 

rms ,nt 

z 

rms , opt 

10 

0.393 

0.149 

0.059 

20 

0.855 

0.323 

0.129 

30 

1.481 

0.560 

0.224 

45 

3.168 

1.198 

0.479 

50 

4.165 

1.574 

0.630 

55 

5.614 

2.122 

0.849 

60 

7.866 

2.973 

1.189 

65 

11.709 

4.426 

1.770 

70 

19.233 

7.270 

2.908 

*This  corresponds  to  the  resonator  mirrors  used 
in  the  interferometric  study.  The  second  column 
gives  the  center-to-edge  values;  the  third  column 
gives  the  rms  OPD  for  an  ultilted  reference 
plane;  the  fourth  column  gives  the  rms  OPD  re- 
lative to  an  optimally-tilted  reference  plane. 


Such  two-dimensional  estimates  would  be  considered  quite  uncertain 
if  standing  alone.  In  later  sections  of  this  report  results  of  more 
detailed  calculations  will  be  presented.  These  rather  simple  sorts  of 
estimates  give  useful  semi-quantitative  guidance  to  the  expected  degree 
of  optical-quality  degradation  for  an  ideally  aligned  system.  On  the 
other  hand,  the  two-dimensional  considerations  give  no  indication  what- 
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TABLE  2.  TWO-DIMENSIONAL  APPROXIMATE  OPD  PREDICTIONS  (IN  MICRONS) 

FOR  A Z CONFIGURATION  TRAVELLING-WAVE  RESONATOR  WITH  R.  = -1481  cm, 

i 

R„  = 3703  cm,  Y = y = 7.5  cm,  FOR  VARIOUS  VALUES  OF  CONVEX- 
2 max 

MIRROR  TILT  ANGLE  6 * 


(deg) 

z 

cp 

z 

rms ,nt 

z 

rms , opt 

10 

0.013 

0.005 

0.002 

20 

0.028 

0.010 

0.004 

30 

0.048 

0.018 

0.007 

45 

0.100 

0.038 

0.015 

50 

0.130 

0.049 

0.020 

55 

0.174 

0.066 

0.026 

60 

0.240 

0.091 

0.036 

65 

0.352 

0.133 

0.053 

68.6 

0.491 

0.185 

0.074 

*This  corresponds  to  a typical  MADS  chemical- 
laser  resonator.  The  second  column  gives  the 
center-to-edge  values;  the  third  column  gives 
the  rms  OPD  for  an  untilted  reference  plane; 
the  fourth  column  gives  the  rms  OPD  relative 
to  an  optimally-tilted  reference  plane. 


ever  of  the  sensitivity  of  optical  quality  to  various  systems  para- 
meters, including  increments  to  mirror  separation  L,  and  concave-mirror- 
tilt-angle  @2  for  a specified  value  of  0^. 


The  two-dimensional  results  give  no  information  about  expected 
shapes  of  interf erograms . Using  information  gathered  from  experiments 
and  three-dimensional  calculations,  a brief  discussion  of  the  behavior 


3 

of  reference-plane  tilt  can  be  given.  In  Figure  2,  the  C y depen- 
dence is  plotted  as  a function  of  z in  the  solid  curve.  If  the  refer- 
ence-plane is  chosen  to  be  that  indicated  by  the  dashed  line,  instead 
of  the  xy  plane  itself,  the  maximum  OPD  will  clearly  be  reduced.  If  it 

2 

is  assumed  that  Z configurations  have  an  x y term  in  the  OPD  series- 
expansion  expression  with  the  same  sign  as  that  of  the  y term,  a 
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TABLE  3.  TWO-DIMENSIONAL  APPROXIMATE  PREDICTIONS  OF  (IN  MICRONS) 

FOR  A Z CONFIGURATION  TRAVELLING-WAVE  RESONATOR  WITH  R1  = -1180  cm, 

R„  = 2640  cm,  Y = y = 10  cm,  FOR  VARIOUS  VALUES  OF  CONVEX- 


R.  = 2640  cm,  Y = y 
2 max 

MIRROR  TILT  ANGLE  Bj* 


(deg) 

z 

cp 

z 

rms ,nt 

z 

rms , opt 

10 

3.057 

0.022 

0.009 

20 

0.124 

0.047 

0.019 

30 

0.216 

0.082 

0.033 

45 

0.468 

0.177 

0.071 

50 

0.619 

0.234 

0.094 

55 

0.839 

0.317 

0.127 

60 

1.184 

0.448 

0.179 

65 

1.778 

0.672 

0.269 

70 

2.952 

1.116 

0.446 

*This  corresponds  to  a typical  CO2  E-beam  laser. 

The  second  column  gives  the  center-to-edge 
values;  the  third  column  gives  the  rms  OPD  for 
an  untilted  reference  plane;  the  fourth  column 
gives  the  rms  OPD  relative  to  an  optimally- 
tilted  reference  plane. 


qualitative  property  of  an  interferogram  of  an  aligned  system  can  be 
predicted  with  the  reference  beam  tilted  as  shown  by  the  dashed  line. 

2 

The  x y term  will  cause  the  OPD  to  increase  for  positive  y when  x is 
varied  in  either  direction  from  zero,  i.e.,  the  surface  for  positive 
y will  curve  upward  with  respect  to  variations  in  x.  Hence,  tilted 
planes  parallel  to  that  indicated  by  the  dashed  line,  but  corresponding 
to  smaller  values  of  z,  will  intersect  the  surface  in  closed  curves. 
Similarly,  for  negative  values  of  y,  tilted  planes  parallel  to  that 
indicated  by  the  dashed  line,  but  corresponding  to  larger  values  of  z, 
will  intersect  the  surface  in  closed  curves.  This  is  the  qualitative 
nature  of  the  interf erograms  for  a suitable  aligned  system  with  tilted 
reference  beam  for  Z configurations. 
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B.  Three-Dimensional  Calculations 

The  major  portion  of  the  theoretical-computational  part  of 
this  study  has  been  concerned  with  developing  and  applying  three- 
dimensional  ray-tracing  computational  methods  and  with  correlating  the 
results  (primarily  in  the  form  of  quasi-interf erograms)  with  the  exper- 
imental inter f erograms . The  development  of  computational  methods  pro- 
gressed from  step-by-step  hand  calculations  to  the  use  of  a programmable 
pocket  calculator,  then  to  a Fortran  program  which  initially  treated 
each  of  the  spherical  mirrors  by  a separate  set  of  equations,  and  finally 
to  a Fortran  program  which  treats  mirror  reflections,  ray  intersections, 
and  coordinate  rotations  by  subroutines  in  such  a way  that  possible 
future  generalizations  to  more  than  two  spherical  mirrors  should  be 
simplified . 

Equations  employed  in  ray  tracing  are  derived  in  Appendix  A.  Their 
present  application  assumes  a collimated  beam  incident  on  the  convex 
mirror.  Given  x and  y coordinates  of  an  incident  ray  measured  relative 
to  the  central  ray,  the  computer  program  calculates  the  path  of  the  ray 
through  the  system  and  determines  (a)  direction  cosines  after  reflection 
from  the  concave  mirror,  measured  in  a coordinate  system  with  the  posi- 
tive z-axis  parallel  to  the  central  ray  and  (b)  the  OPD  between  the  ray 
considered  and  the  central  ray.  The  OPDs  are  of  primary  interest,  and 
are  calculated  for  a position  close  to  the  concave  mirror.  Double- 
precision arithmetic  was  used  to  obtain  reliable  values  of  OPDs. 

The  input  data  include  half-widths  of  the  desired  output  beam  and 
the  number  of  mesh  points  along  x and  y (usually  9 x 9)  in  the  x and  y 
directions.  Coordinates  x and  y of  incident  rays  are  obtained  by 
dividing  corresponding  desired  coordinates  of  output  rays  by  magnifica- 
tions and  , which  are  computed  by  tiie  program  from  input  values 

of  R^,  R2,  and  0^.  Unless  otherwise  specified,  the  program  uses  the 

values  of  mirror  spacing  L and  concave-mirror  tilt  angle  0 which  result 

* 

in  a collimated  output  beam.  Values  of  these  parameters  are  determined 
by  an  initial  portion  of  the  program.  The  first  page  of  computer  out- 
put is  a table  of  calculated  OPDs  over  the  specified  mesh;  the  OPDs  are 
given  in  microns;  mirror  spacing  and  radii  of  curvature  are  given  in 
centimeters. 

As  might  be  expected,  there  is  a good  deal  of  order  in  the  cal- 
culated OPDs;  their  functional  dependence  on  x and  y can  be  well  repre- 
sented by  a relatively  small  number  of  terms  of  a Taylor's  expansion. 

The  number  of  terms  which  are  retained  is  somewhat  arbitrary  and  has 

3 2 

been  increased  from  time  to  time.  Terms  y and  x y are  of  pri- 

03  21  7 

mary  importance.  Values  of  the  coefficients  and  are  not  the 
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same  (as  they  may  be  in  some  contexts  in  which  the  terminology  of  "coma" 

2 2 

is  employed) . The  quadratic  terms  C£q  x and  y can  be  important 

if  the  system  is  not  properly  adjusted.  These  terms  have  been  included 
in  the  function-fitting  process  to  be  able  to  investigate  the  effects 
of  adjusting  L and  e_  to  values  slightly  different  from  their  ideal 

/ 2 2 

values.  A cross-term  C _ x y and  fourth-order  terms  of  the  form 

4 4 ^ 

C/r.  x and  C..  y are  also  considered. 

40  04 

The  computer  program  evaluates  the  series  coefficients  from  the 
table  of  calculated  OPDs  and  prints  these.  It  also  calculates  and 
prints  a table  of  residuals  showing  the  difference  between  the  directly 
calculated  OPDs  and  values  obtained  from  the  truncated  series  expansion; 
these  typically  seem  to  be  of  the  order  of  less  than  one  percent  to  few 
percent  of  the  OPDs  themselves.  The  coefficient  C corresponding  to 

a reference  plane  defined  by  z = C y,  which  minimizes  the  rms  OPD 

is  also  determined  and  printed.  A table  is  printed  (as  the  last  page 
of  output  for  each  case)  of  the  rms  values  of  OPD  for  various  values  of 
reference-beam  tilt  measured  in  units  of  the  optimal  value.  Inspection 
of  this  table  indicates  substantial  reductions  in  the  rms  OPD,  as  would 
be  expected  from  the  discussion  in  connection  with  two-dimensional  cal- 
culations. The  even  more  pronounced  reduction  as  compared  to  peak 
values  of  OPD  can  be  seen  by  comparison  to  those  listed  in  the  OPD 
table. 


If  L and/or  6^  are  set  to  values  incremented  slightly  from  their 

proper  values,  there  will  be  curvature  in  either  or  both  the  x and  y 
directions,  i.e.,  C^q  and/or  will  differ  appreciably  from  zero. 

Expected  values  of  these  coefficients  can  be  derived  analytically  and 
compared  with  results  of  the  function-fitting  process  applied  to  the 
ray-tracing  OPD  calculations.  This  is  discussed  in  Appendix  B.  Agree- 
ment between  analytically  predicted  and  "function-fitting"  calculated 
values  is  felt  to  be  a confirmation  of  proper  functioning  of  the  computer 
program. 


The  primary  experimental  data  regarding  optical  quality  from  the  >\ 

tilted  spherical  mirror  resonator  are  in  the  form  of  interferograms. 

While  providing  very  precise  and  detailed  information  about  the  optical 
systems,  interferograms  do  not  immediately  (without  an  intervening 
measurement  and  data-reduction  process)  lead  to  quantities  predicted  by 
the  calculations.  Therefore,  it  seemed  useful  to  convert  the  calculated 
results  into  the  form  of  computer-simulated  quasi- interferograms . For 
this  purpose  the  computer  generates  a fairly  large  array  (usually 
49  x 49)  of  OPDs  obtained  by  evaluation  of  the  truncated  series  expansion 
at  each  of  the  array  points.  The  array  of  numerical  values  is  converted 
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to  an  array  of  alpnabetic  characters  or  blanks  and  used  to  generate  a 
one-page  printer  plot  whose  general  apperance  simulates  that  of  an 
interf erogram.  Successive  characters  of  the  alphabet  correspond  to 
incremental  OPDs  of  one  wavelength  (0.6328  pm).  The  interspersing  of 
blanks  with  letters  of  the  alphabet  has  the  result  that  (at  least  in 
the  central  portion  of  the  plot)  light  areas  are  interspersed  with  dark 
areas  and  properly  simulate  an  interferogram.  With  only  26  characters 
available,  the  resolution  is  limited.  As  a result,  it  often  happens 
that  casual  inspection  of  outer  portions  of  such  plots  suggests  spurious 
fringe  pattern,  because  the  eye  tends  to  interpret  contiguous  light 
areas  as  all  belonging  to  a light  fringe.  This  difficulty  can  be 
avoided  by  concentrating  on  dark  areas  and  keeping  in  mind  that  a par- 
ticular dark  fringe  is  always  represented  by  the  same  letter  of  the 
alphabet.  In  one  sense,  the  quasi-interferograms  give  more  information 
than  genuine  interferograms  because  the  sign  and  the  magnitude  of  the 
OPD  relative  to  the  center  are  represented. 

The  computer  program  will  optionally  produce  multiple  quasi- 
interferograms  corresponding  to  a specified  set  of  reference-beam  tilts 
from  a single  set  of  OPD  calculations  (Figures  3 through  14) . The 
observed  progression  of  shapes  as  tilt  is  systematically  varied  has 
proven  interesting  and  useful  because  reference-beam  tilt  can  also  be 
varied  experimentally.  Inspection  of  quasi-interferograms  leads  to  some 
interesting  observations  for  a properly-aligned  system  (which  was  treated 
computationally  prior  to  investigating  effects  of  incremented  values  of 
L and  0^) : 

1)  There  is  symmetry  in  x,  as  would  be  expected  since  the  entire 
optical  system  is  assumed  to  have  such  symmetry. 

2)  For  Z configurations  and  small-to-moderate  values  of  0 , the 

pattern  is  symmetric  in  y to  a fairly  high  degree  of  approximation  (the 
OPDs  are  approximately  antisymmetric) . 

3)  For  Z configurations  there  is  a regular  progression  of  shapes 
as  reference-beam  tilt  is  increased  from  zero  to  its  optimal  value  and 
beyond.  Specifically,  the  pattern  changes  from  a single  central  oval 

to  a pair  of  ovals  (symmetrically  spaced  about  y = 0)  which  move  further 
apart  and  are  separated  by  an  increasing  number  of  fringes. 

These  observed  properties  are  useful  in  experimental  adjustment  of 
L and  0 2»  because  the  symmetry  in  y is  a rather  good  "signature".  Cal- 
culations made  with  values  of  L and/or  02  which  are  incremented  from 

their  ideal  values  produce  quasi-interferograms  which  are  noticeably 
lacking  in  the  y-symmetry  property. 

From  one  point  of  view,  the  tilted-spherical-mirror  resonator 
concept  can  be  described  as  a scheme  in  which  a large  degree  of  astig- 
matism is  introduced  by  the  tilt  of  the  convex  mirror,  and  is  cancelled 
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by  suitable  choice  of  mirror  spacing  and  the  tilt  of  the  concave  mirror. 
Quasi- inter ferograms  can  be  used  to  illustrate  cancellation  of  curva- 
ture in  either  the  x or  y coordinates  (simultaneous  cancellation  along 
both  coordinates  is  achieved  by  the  properly  adjusted  tilted-mirror 
resonator).  For  this  purpose,  a small  value  (10°)  was  chosen  for  0^  so 

that  the  effects  to  be  illustrated  would  not  be  masked  by  other  aberra- 
tions. The  sequence  of  quasi- interferograms  is  presented  in  Figures  15 
through  19.  The  basic  case,  i.e.,  the  properly  adjusted  system,  shows 
very  little  aberration  of  any  sort;  hence,  it  produces  a rather  struc- 
tureless quasi-interferogram.  When  a small  increment  AL  (labelled 
XINCL)  in  mirror  spacing  is  introduced,  curvature  of  the  same  sign  is 
introduced  along  both  x and  y axes;  the  characteristic  hill  or  valley 
shape  is  evident  in  Figure  16.  Instead,  if  L is  kept  at  its  proper  value 
but  a small  increment  A0^  (labelled  XINCT)  to  concave-mirror  tilt  is 

introduced,  curvatures  with  opposite  sign  along  the  x and  y axes  are 
introduced.  This  leads  to  the  characteristic  "X"  pattern  representing 
a saddle-point  type  of  surface,  as  seen  in  Figure  17.  For  a specified 
value  of  AL  (or  of  A@2)  an  increment  of  the  other  parameter  can  be 

chosen,  such  as  to  cancel  the  curvature  in  either  the  x or  y dimension. 
Such  cases  are  presented  in  the  Figure  18  and  19  plots,  which  show 
vertical  and  horizontal  fringes.  It  should  be  emphasized  that  such 
cancellation  could  be  obtained  if  there  were  additional  sources  of 
astigmatism  present  in  the  system.  When  one  is  free  to  adjust  L and 

curvature  can  be  eliminated  in  the  x and  y dimensions;  i.e.,  coeffi- 
cients C2Q  and  Cq2  are  reduced  to  zero.  This  is  the  criterion  which 

defines  the  tilted-spherical-mirror  resonator  concept.  It  should  empha- 
size that  this  cancellation  can  also  be  obtained  if  other  sources  of  astig- 
matism are  present.  This  small  value  of  0^  was  deliberately  chosen  to 

illustrate  avoiding  complicating  the  plots  with  substantial  amounts 

2 2 

of  aberrations  other  than  those  associated  with  x and  y expansion 
terms.  For  more  realistic  cases  with  larger  values  of  0 , the  quasi- 

interf erograms  do  not  have  such  a simple  appearance. 

A number  of  calculations  were  made  for  0^  = 45°  for  various  non- 
zero values  of  AL  and  ^ (as  well  as  reference-beam  tilts);  results 

are  shown  in  Figures  20  through  43.  The  coefficient  of  either  the 
2 2 

x or  y term  can  also  be  reduced  to  zero  by  suitable  choice  of  one  of 
the  increments  when  the  other  is  specified.  Reduction  to  (approxi- 

2 

mately)  zero  of  the  x term  shows  up  most  clearly;  the  plot  shows 

that  there  is  very  little  OPD  variation  along  x for  y = 0.  When  the 
2 

coefficient  of  y is  reduced  to  zero,  the  corresponding  quasi 
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interferogram  is  somewhat  more  involved  because  there  is  still  an  appre- 
ciable amount  of  aberration  present  associated  with  other  expansion  terms. 
What  is  especially  noticeable  is  that  the  y-symmetry  is  decidedly 
absent;  this  emphasizes  usefulness  of  the  y-symmetry  property  in  making 
experimental  adjustments.  If  the  signs  of  AL  and  A0„  are  reversed,  the 

2 1 
y coefficient  still  vanishes,  there  is  still  a lack  of  y-symmetry,  and 

the  plot  is  essentially  just  the  mirror  image  of  the  plot  obtained 

without  reversal  of  signs.  This  mirror-image  property  applies  for  any 

pair  of  magnitude  of  AL  and  In  the  special  case  of  zero  magnitude 

of  these  increments,  the  mirror  image  of  the  pattern  is  simply  (to  a 
rather  good  approximation)  identical  with  the  pattern  itself;  i.e.,  the 
y-symmetry  property  is  characteristic  of  proper  adjustment  of  L and  62* 

The  y-symmetry  property  is  not  exact,  though  it  seems  to  hold  well 
for  Z configurations  (which  are  of  greatest  interest  because  they  have 
fewer  aberrations  than  U configurations)  up  to  large  values  of  0^. 

Presumably,  this  property  results  from  the  fact  that  the  principal 
terms  in  the  series  expansion  of  the  OPD  are  of  odd  power  in  y 

3 2 

(specifically,  they  are  the  y and  x y terms). 

For  cases  where  comparisons  have  been  made,  the  three-dimensional 
ray-tracing  calculations  predict  optical-quality  parameters  which  are 
not  greatly  different  from  those  obtained  by  the  much  simpler  approach 
based  on  the  two-dimensional  approximation.  Some  specifics  will  be 
given  in  the  following  paragraphs. 

For  the  representative  MADS-type  chemical-laser  resonator  with 
0^  = 68.6°  and  output-beam  dimensions  of  3 * 15  cm,  the  detailed  cal- 
culations give  a center-to-peak  OPD  along  the  y-dimension  (for  x = 0) 
of  approximately  0.49  pm,  in  agreement  with  the  two-dimensional  pre- 
diction of  Table  2.  The  rms  value  of  the  OPD  for  an  untilted  reference 
plane  is  calculated  to  be  0.188  pm  as  compared  to  a two-dimensional 
prediction  of  0.185  pm.  The  minimized  (with  respect  to  reference-plane 
tilt)  rms  value  of  OPD  is  calculated  to  be  0.075  pm  as  compared  to  the 
two-dimensional  prediction  of  0.074  pm.  The  agreement  is  rather  good. 
This  is  partly  due  to  the  fact  that  the  flow  dimension  is  small  com- 
pared to  the  cross-flow  dimension,  so  that  the  resonator  is,  in  fact, 
something  of  a strip. 

For  the  typical  C02  resonator  at  9^  = 55°,  detailed  calculations 

give  a center-to-peak  OPD  along  the  y-dimension  (for  x = 0)  of  approxi- 
mately 0.837  pm  and  a center-to-corner  OPD  of  approximately  1.49  pm; 
by  comparison  the  two-dimensional  center-to-edge  prediction  of  Table  3 
is  0.839  pm.  The  rms  value  of  OPD  for  an  untilted  reference  plane  is 
0.448  pm  compared  to  a two-dimensional  prediction  of  0.317  pm.  The 
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minimized  rms  OPD  is  calculated  to  be  0.169  um  as  compared  to  a two- 
dimensional  prediction  of  0.127  pm. 


For  the  mirror  parameters  associated  with  the  experimental  studies, 
calculations  have  been  carried  out  for  a range  of  values  of  0^,  for 

which  the  results  are  listed  in  Table  4. 


TABLE  4.  THREE-DIMENSIONAL  RAY-TRACING  PREDICTIONS  (IN  MICRONS) 

OF  OPD  FOR  A Z CONFIGURATION  TRAVELLING-WAVE  RESONATOR  WITH 

R,  = -290  cm,  R„  = 675  cm,  x = y =7.62  cm,  FOR  VARIOUS 
1 2 max  max 

VALUES  OF  CONVEX-MIRROR  TILT  ANGLE  Bj* 


6X  (deg) 

z 

cp 

z 

rms,nt 

z 

rms, opt 

10 

0.393 

0.230 

0.091 

20 

0.855 

0.496 

0.194 

30 

1.481 

0.848 

0.330 

45 

3.171 

1.764 

0.678 

50 

4.171 

2.287 

0.876 

55 

5.626 

3.031 

1.157 

60 

7.892 

4.159 

1.582 

70 

19.450 

9.529 

3.667 

*This  corresponds  to  the  resonator  mirrors  used 
in  the  interferometric  study.  The  second 
column  gives  the  center-to-edge  values  (taken 
along  y for  x = 0) ; center-to-corner  values 
are  somewhat  larger  but  are  not  listed  here. 
The  third  column  gives  the  rms  OPD  for  an 
untilted  reference  plane;  the  fourth  column 
gives  the  rms  OPD  relative  to  an  optimally- 
tilted  reference  plane. 


These  results  may  be  directly  compared  with  the  two-dimensional  predic- 
tions for  the  same  case  as  given  in  Table  1.  The  three-dimensional 
predictions  of  the  optimized  rms  OPD  range  from  only  25%  to  50%  larger 
than  the  two-dimensional  predictions,  even  for  this  case  which  is  a 
resonator  with  larger  aberrations. 
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The  adjustment  of  L and  when  R^  and  are  fixed  and  0^  is  set, 

is  a central  feature  of  this  resonator  concept.  As  alluded  to  earlier, 
a significant  aspect  of  the  recent  studies  has  been  development  of  cri- 
teria and  procedures  for  making  the  appropriate  adjustments.  The  question 
naturally  arises  as  to  the  quantitative  degree  of  sensitivity  of  rms 
OPD  to  variations  in  L and  62*  Some  numerical  results  are  plotted 

in  Figures  44  through  46.  It  was  found  that  the  numerically-determined 
values  are  in  fairly  good  agreement  with  predictions  that  can  be  made  on 
the  basis  of  the  intially  somewhat  questionable  assumption  that  aberra- 
tions othe**  than  those  of  quadratic  type  are  insensitive  to  adjustment 
of  L and  62 , and  that  these  other  aberrations  are  odd  in  y and  hence 

"noninteracting"  with  the  quadratic  aberrations  introduced  by  variations 
in  L and/or  0„.  Following  the  approach  of  Appendix  B,  the  ratio  of  the 

coefficients  C?^  and  C ^ to  AL  and/or  can  then  be  determined.  With 

these  assumptions,  values  of  rms  OPD  can  be  predicted  as  a function  of 
AL  for  A02  = 0 and  as  a function  of  A@2  for  AL  = 0.  These  predictions 

are  plotted  as  the  smooth  curves  in  Figures  45  and  46.  The  numerically- 
determined  values  agree  satisfactorily  with  the  curves.  The  OPD  con- 
tributions from  increments  to  L and  to  @2  are  not  additive  (though  the 

magnitude  of  the  combined  effect  of  specified  increments  to  both  para- 
meters can  be  predicted) . 

For  the  system  being  studied  experimentally  at  moderate  values  of 
3^  (45°  and  55°),  it  appears  that  L should  be  adjusted  to  within  a few 

millimeters,  and  ©2  to  within  a few  milliradians  in  order  that  wave- 

curvature  effects  not  increase  the  overall  aberration  substantially 
beyond  the  minimal  value  which  is  inherent  in  the  tilted-spher ical  mirror 
design.  These  are  not  very  stringent  restrictions  on  experimental 
adj  ustments . 

III.  EXPERIMENTAL  RESULTS 

The  results  of  this  section  agree  quantitatively  with  the  computa- 
tional results  obtained  in  the  three-dimensional  calculations  of  wave- 
front  OPDs.  A set  of  mirrors  was  chosen  to  give  a moderate  number  of 
fringes  across  the  output:  aperture  and  could  be  compared  directly  to 
calculated  quasi-interferograms.  This  test  represented  an  extreme  case 
compared  to  the  aberrations  which  would  be  present  in  the  proposed  high 

aspect  ratio  resonator  mirrors^  but  provided  a direct  correlation 
between  theory  anc'  experiment. 


6Ibid. 
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Two  types  of  resonator  concepts  were  tested,  the  Z configuration 
and  the  U configuration.  The  schematics  of  these  configurations  are 
shown  in  Figures  1(a)  and  1(b).  Most  of  the  interferograms  were  taken 
with  the  Z configuration  experimental  setup.  Only  two  photographs  were 
taken  in  the  U configuration.  All  agreed  well  with  the  analytical 
predictions  of  OPD. 

Schematics  of  the  experimental  setup  for  the  Z and  U configurations 
are  shown  in  Figures  47(a)  and  47(b).  In  both  cases  a plane  wavefront 
is  produced  by  passing  a helium-neon  laser  beam  through  a spatial  filter 
(microscope  objective  plus  pinhole)  and  collimating  the  resultant  spher- 
ically divergent  wavefront  with  a high  quality  objective  lens  of  6-in. 
diameter.  This  plane  wavefront  is  then  split  into  two  components  with 
a glass  beamsplitter.  One  component  (reference  beam)  is  reflected 
directly  into  a focusing  lens  to  a camera.  The  other  portion  (test 
beam)  passes  through  the  beamsplitter  to  the  convex  resonator  mirror, 
reflects  to  the  concave  resonator  mirror,  and  is  directed  back  through 
the  beamsplitter  into  the  camera  by  means  of  appropriately  positioned 
mirror  flats. 

To  obtain  the  best  quality  interferograms,  an  alignment  procedure 

for  adjustment  of  the  resonator  mirrors  is  necessary.  First,  the  angles 

0. , 0 and  the  separation  between  the  mijprors  is  setup  in  gross  agree- 

1 ^ 1 
ment  with  the  theoretical  predictions  (Figure  1) . With  the  reference 

beam  blocked,  the  beam  passing  through  the  resonator  and  focusing  lens  is 

observed.  The  concave  resonator  mirror  is  tilted  (variation  of  0 ) with 

simultaneous  compensation  of  the  turning  flat  mirror  to  keep  the  beam 
on  the  optical  axis  of  the  focusing  lens.  The  optimal  setting  of  is 

obtained  with  the  sharpest  focal  point. 

For  any  given  value  of  0^  there  is  only  one  value  of  6^  and  L which 

minimize  the  aberrations  introduced  in  the  tilted  spherical  mirror  con- 
figurations. The  mirror  separation  L is  not  as  critical  as  0^  to  set 

up  experimentally.  No  variation  from  the  theoretically  determined  value 
was  necessary. 

The  interferograms  were  obtained  by  superimposing  the  focal  points 
of  the  two  beams  from  the  beamsplitter  (by  tilting  the  turning  flat). 
There  is  a certain  arbitrariness  in  this  adjustment—  the  interferogram 
obtained  depends  strongly  upon  the  angle  between  the  reference  beam  and 
the  test  beam  (wavefront  tilt).  This  variability  made  it  difficult  to 
compare  exactly  the  quasi- interferograms  obtained  analytically  with  the 
experimental  interferograms. 

7Ibid. 
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Figures  48  through  54  illustrate  the  effect  of  variation  of  wave- 
front  tilt  (reference  beam  tilt)  upon  the  appearance  of  the  interfero- 
grams.  A direct  side-by-side  comparison  of  theory  and  experiment  is 
given  in  many  of  these  and  the  following  cases  shown.  The  agreement  is 
seen  to  be  very  good.  These  cases  are  for  nominal  alignment  with 
0^  = 45°  in  the  Z configuration.  A vignetting  of  the  6-in.  square 

theoretical  OPD  plots  is  a result  of  the  concave  mirror  being  tilted 
with  respect  to  the  optic  axis  (Figure  1) . The  mirror  is  round  and  when 
viewed  from  the  output  end  (2b),  it  produces  a test  beam  having  an 
elliptical  cross  section  of  6-in.  major  axis  by  6 cos  in*  minor  axis. 

This  has  been  taken  into  account  by  not  drawing  an  interferogram  in  the 
portions  of  the  theoretical  plots  not  seen. 

Effects  of  misalignments  (ALjAS^)  for  0^  = 45°  are  shown  in 

Figures  55  through  58.  Figure  55  shows  the  result  of  a change  in  sepa- 
ration (AL)  between  the  convex  and  concave  mirrors  (closer  together  by 
1 cm)  . Figures  56  and  57  result  from  varying  approximately  ±7  mrad 

(A02)  about  the  correct  angular  position.  The  nirror  symmetry  intro- 
duced about  the  vertical  axis  for  A62  of  +7  mrad  is  evident  in  the  inter- 

ferograms.  An  increase  in  deviation  of  ©2  of  14  mrad  increases  the 
aberrations  as  seen  by  the  increased  number  of  fringes  in  Figure  58. 

The  U configuration  was  tested  for  0^  = 45°.  Figures  59  and  60 

are  interferograms  obtained  for  two  settings  of  reference  beam  tilt. 

The  aberrations  increase  with  increasing  0^  as  shown  in  Figures  61-74 

for  0^  = 70°.  Figures  61-64  show  a variation  in  reference  beam  tilt 

for  the  nominally  aligned  position  of  resonator  mirrors.  A comparison 
of  this  set  of  figures  with  the  similar  set  at  0 = 45°  (Figures  48 

through  54  show  that  the  essential  features  of  the  interferograms  are 
the  same,  but  that  the  aberrations  are  worse  (more  fringes)  in  the  70° 
case . 

Variations  in  AL  and  A©2  for  6^  = 70°  show  the  same  general 

features  as  for  A0^  = 45°.  Figures  65  and  66  are  interferograms  taken 

for  AL  = ±2  cm  variation  about  the  nominal  separation  in  spherical 
mirrors.  Figures  67  and  68  show  = ±7  mrad.  Figures  69  and  70  give 

a variation  of  = ±14  mrad.  The  mirror  symmetry  of  these  interfero- 
grams in  AL  and  kQ ^ is  evident  in  the  figures.  These  interferograms 
were  taken  in  the  Z configuration. 
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The  interf erograms  for  the  U configuration  show  worse  phase  front 
errors  than  the  Z configuration.  Figures  71  and  72  are  interf erograms 
at  0^  = 70°  for  two  cases  of  reference  beam  tilt.  These  interferograms 

as  well  as  the  theoretical  counterparts  are  not  symmetric  about  the 
vertical  axis.  These  are  due  to  the  higher  order  even  terms  contrib- 
uting to  the  phase  error  in  the  series  expansion  of  the  OPD. 

Figures  73  and  74  are  a return  to  the  0^  = 70°,  Z configuration  and 

show  a tilt  of  the  concave  mirror  out  of  the  plane  of  rotation  of 
0^  and  6^  (in  the  vertical  direction  of  +7  mrad) . The  mirror  symmetry 

of  the  interferograms  about  the  horizontal  axis  for  these  misalignments 
is  seen. 


IV.  SUMMARY  AND  CONCLUSIONS 

Advantages  of  this  relatively  new  type  of  resonator  have  been, 
noted  earlier;  these  include  asymmetric  magnification  (M  ^ M ) com- 
bined with  collimated  output,  a considerable  measure  of  experimental 
and  design  flexibility,  and  the  relative  ease  of  fabrication  and  adjust- 
ment of  spherical-figure  mirrors  (as  used  in  the  present  designs)  as 
constrasted  to  the  case  of  toroidal  mirrors  (which  are  otherwise  used  to 
achieved  asymmetric  magnif ication) , The  disadvantages  of  tilted- 
spherical-mirror  resonators  is  the  additional  deterioration  of  output 
optical  quality  resulting  from  the  substantial  mirror  tilt  angles  which 
are  involved;  a quantitative  investigation  of  this  effect  seemed  clearly 
in  order.  Investigation  of  the  optical  quality  is  proceeding  in  stages. 
A preliminary  interferometric  investigation  was  already  reported  in 

8 9 

connection  with  the  initial  examination  of  the  resonator  concept.  ’ 

This  study  has  considered  in  a rather  detailed  manner  the  effect  on 
optical  quality  resulting  from  one  complete  passage  from  convex  to 
concave  mirror  (essentially  a conventionally  aligned  travelling-wave 
configuration)  of  an  input  collimated  beam.  Comparison  of  computational 
and  experimental  results  has  been  facilitated  by  presenting  computed 
results  in  a format  similar  in  appearance  to  interferograms.  Computa- 
tional and  experimental  results,  including  shapes  of  interferograms  and 
their  behavior  as  a function  of  reference-beam  tilt  angle,  incremental 
mirror  separation,  and  incremental  concave  mirror  tilt  angle,  are  in 
satisfactory  agreement.  The  considerable  reduction  in  root-mean-square 
OPD  which  results  from  using  a slightly  tilted  reference  plane  has  been 
confirmed;  this  simply  implies  that  in  a practical  device,  a major 


8Ibld. 
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portion  of  the  distinctive  tilted-mirror  effect  consists  of  a slight  angu- 
lar displacement  of  far-field  peak  intensity  relative  to  that  of  the  nomi- 
nal position.  Small  angular  pointing  errors  are  correctable  by  a tracker. 

The  experimental  and  computational  studies  have  both  clearly  shown 
that  the  Z configuration  is  to  be  preferred  over  the  U configuration  as 
regards  output  optical  quality.  The  principal  tilted-spherical-mirror 
aberrations  in  a properly  adjusted  system  a^ise  from  contributions  to 

3 2 

optical  path  differences  which  are  proportional  to  y and  to  x y.  For 

3 

a high-aspect-ratio  resonator,  only. the  y contribution  will  be  of 
substantial  significance;  the  magnitude  of  this  contribution  was  found 
3 2 

to  scale  as  y /R  . Formulas  are  given  for  determining  the  magnitude  of 
this  contribution  in  a one-transverse-dimension  approximation  without 
extensive  numerical  computations,  though  three-dimensional  ray-tracing 
calculations  are  required  to  evaluate  the  detailed  effects.  Tables  are 
given  of  results  of  approximate  and  detailed  calculations  for  selected 
systems . 

Certain  additional  aspects  of  this  general  class  of  resonators 
could  be  investigated  by  extensions  of  methods  used  in  the  present 
study.  First,  the  effect  of  repeated  propagation  of  a beam  through  the 
resonator  (satisfying  the  round-trip  consistency  condition)  should  be 
investigated;  this  will  include  effects  of  decentered  alignment. 

Secondly,  and  more  importantly,  the  studies  should  be  extended  to  some- 
what more  involved  configurations  than  those  considered  here.  In  par- 
ticular, an  investigation  needs  to  be  made  of  standing-wave  tilted- 
spherical-mirror  configurations,  which  involve  a more  complex  optical 
path  including  a turning  flat  or  inverter  associated  with  each  of  the 
two  spherical  mirrors.  In  such  an  investigation  an  effort  should 
naturally  be  made  to  find  methods  of  optimizing  optical  quality  as  well 
as  achieving  other  desirable  features  (e.g.,  reduced  sensitivity  to 
laser-medium  inhomogeneities) . Distances  between  turning  flats  or 
inverters  and  their  associated  spherical  mirrors  are  additional  experi- 
mental variables,  suitable  adjustment  of  which  may  possibly  lead  to 
improvement  of  optical  quality  as  compared  to  the  travelling-wave  con- 
figuration, as  well  as  providing  variations  in  magnifications  M and  M . 

x y 

An  especially  desirable  goal  in  this  connection  is  finding  combinations 
of  parameters  such  that  aberration  effects  are  reduced  to  very  small 
magnitudes.  While  there  is  no  compelling  reason  for  expecting  this  to 
be  possible,  it  cannot  be  ruled  out  in  advance,  particularly  because 

3 

in  typical  chemical- laser  configurations,  only  the  y aberration  term 
should  be  of  substantial  significance. 

The  effects  of  beam  inverters  associated  with  either  or  both  spher- 
ical mirrors  in  a standing-wave  device  should  be  investigated,  not  only 
as  a possible  means  of  reducing  empty-resonator  aberrations  but  also 
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because  beam  inversion  could  be  quite  helpful  in  reducing  deleterious 
effects  of  spatial  nonuniformity  of  gain  and  index  of  refraction  in  the 
laser  gain  medium.  It  is  also  desirable  to  investigate  the  latter 
effects  separately,  perhaps  by  wave-optics  calculations  and  by  experi- 
ments. The  effects  of  beam  inversion  in  travelling-wave  resonators 
should  also  be  considered;  it  should  be  noted  that  inversion  in 
travelling-wave  resonators  can  be  accomplished  with  an  odd  number  of 
flaL  turning  mirrors  without  the  necessity  of  special  inverting  optical 
components. 


Y3  aberration 
TERM 


REFERENCE 

BEAM 


NOTE:  THE  RMS  VALUE  OVER 
THE  APERTURE  (RANGING  FROM 
-Y  to  +Y)  CAN  BE  MINIMIZED  BY 
REFERENCE  BEAM  TILT  (RBT). 


Figure  2.  Dependence  of  the  y aberration  term. 
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Figure  3.  Computer-generated  quasi-inter ferogram  Z configuration, 
= 45°,  nominal  values  of  L and  e^,  reference-beam  tilt  (rbt) 

(optimal  value  of  RBT)  - 0.222  (-0.370),  in  units  of  10-4  rad. 
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Figure  4.  Z configuration,  0^  = 45°,  RBT  = 0.148  (-0.370), 
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Figure  5.  Z configuration,  01  = 45°,  RBT  = 0.074  (-0.370), 
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Figure  6.  Z configuration,  0^  = 45°,  RBT  = 0.000  (-0.370), 
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Figure  7.  Z configuration,  0.  = 45°,  RBT  = -0.074  (-0.370), 
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Figure  8.  Z configuration,  9^^  = 45°,  RBT  = -0.148  (-0.370). 
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Figure  9.  Z configuration,  0X  = 45°,  RBT  = -0.222  (-0.370). 
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Figure  11.  Z configuration,  0^  = 45°,  RBT  = -0.370  (-0.370). 
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Figure  12.  Z configuration,  0- 
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-0.444  (-0.370). 
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Figure  13.  Z configuration,  0^  = 45°,  RBT  = -0.518  (-0.370) . 
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Figure  14, 
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Figure  15.  Z configuration,  0^  = 10°,  RBT  = 3,  nominal  values 
of  L and  > i»e«>  AL  = 0,  = 0- 
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Figure  16.  Z configurations,  0^ 
A02  = 0. 


10°,  RBT  = 0,  AL  = 1.3097  cm, 
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Figure  17.  Z configuration,  6^  = 10°,  RBT 


0,  AL  = 0, 


A62  = -0.04233  rad. 
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Figure  18.  Z configuration,  0 = 10°,  RBT  = 0,  AL  = 1.3097  cm, 

A©2  = -0.04233  rad.  The  curvature  along  the  y-coordinate  has 
been  reduced  to  approximately  zero. 
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Figure  19.  Z configuration,  0^  = 10°,  RBT  = 0,  AL  = -1.333  cm, 
A0£  = -0.04233  rad.  The  curvature  along  the  x-coordinate  has 
been  reduced  to  approximately  zero. 
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Figure  20.  Z configuration,  0^  = 45°,  RBT  = 0 (-0.379). 
AL  = 1.0  cm,  A0O  = 0. 
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Figure  21.  Z configuration,  0-  = 45c 
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Figure  22.  Z configuration,  0^  = 45°,  RBT  = 0 (-0.361), 


L = -1.0  cm,  0, 


0. 
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Figure  23.  Z configuration,  0^  = 43°,  RBT  = -0.361  (-0.361), 
AL  = -1,0  cm,  A0£  = 0. 
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Figure  24.  Z configuration,  0^  = 45°,  RBT  = 0 (-0.389), 


AL  = 2.0  cm,  A0£  = 0. 
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Figure  25.  Z configuration,  9^  = 45°,  RBT  = -0.389  (-0.389), 


AL  = 2.0  cm,  A0£  = 0. 
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Figure  26.  Z configuration,  0^  = 45°,  RBT 


-0  (-0.351) 


AL  = -2.0  cm,  = 0. 
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Figure  27.  Z configuration,  0.  = 45°,  RBT  = -0.351  (-0.351), 
AL  = -2.0  cm,  A0^  = 0. 
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Figure  28.  Z configuration,  0 = 45°,  RBT  = 0 (-0.392), 


AL  = 0,  A02  = 0.007  rad. 
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Figure  29.  Z configuration,  0 = 45°,  RBT 
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Figure  30.  Z configuration,  0^  = 43°,  RBT  = 0 (-0.348),  AL 


0, 


A02  = -0.007  rad. 
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Figure  31.  Z configuration,  9^  = 45°,  RBT  = -0.348  (-0.348), 
AL  = 0,  A92  = -0.007  rad. 
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Figure  32.  Z configuration  0 


x 45°,  RBT  = 0 (-0.393),  AL  = 1.0  cm, 
The  curvature  along  the  x-coordinate  has  been 
reduced  to  approximately  zero. 


A0^  = 0.00428  rad 


57 


I 


QQQQQ  PPPP 

000 

NNN 

QQQQQQ  PPPP 

000 

NNN 

QQQQ  PPPP 

0000 

NNN 

QQ  PPPPP 

000 

NNN 

PPPPP 

0000 

NNN 

PPPPPP 

0000 

NNNN 

PPPPPP  0000 

NNN 

PPPPPP  00000 

NNNN 

PPP  00000 

NNNN 

000000 

NNNN  M 

0000000 

NNNNN  M 

000000000 

NNNNNN  MM 

0000000  NNNNNN 

MM 

000  NNNNNNN 

MMM 

NNNNNNNN 

MMMM 

NNNNNNNNNN 

MMMMM 

NNNNNNNNNNNN 

MMMMM 

NNNNNNNNNNNN 

MMMMMM 

INNNNNNNNN 

MMMMMMM 

SS  TTT  uuuu  VVVVVVV  VVVVVVV  UUUU  T T T SS 

R SSS  TTT  UUUUU  VVVVVVVVVVV  VVVVVVVVVVV  UUUUU  TTT  SSS  R 

RR  SSS  TTT  UUUUU  V VVVV VVV VVV V VVV VV  UUUUU  TTT  SSS  R R 

RRR  SSS  TTTT  UUUUUUUU  UUUUUUUU  TTTT  SSS  R R R 

RRR  SSS  TTTTT  UUUUUUUUUUUUUUUUUUUUUUUUUUU  TTTTT  SSS  RRR 

RRR  SSSS  TTTTTTT  UUU  TTTTTTT  SSSS  RRR 

QQ  R RR  R SSSSS  TTTTTTTTTTT  TTTTTTTTTTT  SSSSS  RRRR  Q Q 

QQQG  RRRR  SSSSSS  TT TT TT T TT T T TT T T T T T T SSSSSS  RRRR  0000 

00  Q 0 RRRRR  SSSSSSSS  SSSSSSSS  RRRRR  Q QQ  Q 

Q00QQ  RRRRRR  S S S S S S S S S S S S S S S S S S S S S S S S S S S S S RRRRRR  QQQQQ 

PP  QQQQQ  RRRRRRRR  RRRRRRRR  QQQQQ  PP 

ppppp  QQQQQQ  RR RRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRR  QQQQQQ  PPPPP 

PPPPP  QQQQQQQQ  R R R R R R R R R R R R R QQQQQQQQ  PPPPP 

Pppppp  QQQQQQ QQQQQQ  QQQQQQQQQQQQ  PPPPPP 

000  PPPPPPPP  QQQQQQQQQQQQQQQQQQQQQQQQQQQQQ  PPPPPPPP  000 

000000  PPPPPPPPPP  PPPPPPPPPP  000000 

00000000  PPPPP PPPPPPPPPPPPPPPPPPPPPPPPPPPPP PPPPP PPPP  00000000 

0000000000  PPPPPPPPPPPPPPPPPPP  0000000000 

N 000000000000000  000000000000000  N 

NNNN  0000000000000000000000000000000000000000000000000  NNNN 

NNNNNNNNNN  00000000000000000  NNNNNNNNNN 

NN NNNN NN NNNNNNNNNNN  NNNNNNNNNNNNNNNNNNN 

NNNNN NNNN NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN NNNN NNNNNNNNNN 
NNNNNNNNNN NNNN NNNNNNNNN NNNNNNNNNNN NNNNNNNNN NNNNN NNNN NNNN NNNNN 


MMMMMWMMMNMMMWMMMWMWWMMMNMMMWNNWNNW 

NN  NN 

NNNNN  NNNNN 

NNNNNNNNN  WMMM NNNNN NNNN 

NNNNNNNNN  PMMNMMMMMMHMMMNMMMWMPMMMMMMMMMPMM  NNNNNNNNN 

NNNNNNNNN  NNNNNNNNN 

0000  NNNNNNNNNN  NNNNNNNNNN  0000 

000000  NNNNNNNNNN N NNNNNNNNNNN  000000 

PP  000000  NNNNN NNNNNNNNNNNN  NNNNNNNN NNNNN NNNN  000000  PP 

PPPP  0000000  NNNNNN NNNNNNNNNNN NNNNNNNNNN  0000000  PPPP 

QQ  PPPPP  00000000  00000000  PPPPP  QQ 

QQQQ  PPPPP  00000000000000000000000000000000000  PPPPP  QQQQ 

RRRR  QQQQ  PPPPPPP  PPPPPPP  QQQQ  RRRR 

SS  RRRR  QQQQQ  PP P P P P P P P P P P P P P P P P P P P P P P P P P P P P P P P QQQQQ  RRRR  SS 

SSS  RRRRR  QQQQQQQQ  QQQQQQQQ  RRRRR  SSS 

TTT  SSSS  RRRRRR  Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q Q RRRRRR  SSSS  TTT 


UUU 

TTTT 

SSSSS 

RRRRRRRRPRRRRRRRRRRRRRRRRRRRR 

SSSSS 

TTTT 

UUU 

VV  V 

UUU 

TTTTT 

sssssssss 

SSSSSSSSS 

TTTTT 

UUU 

VVV 

MM 

vvv 

UUUU 

TTTTTTT 

TTTTTTT 

uuuu 

vvv 

ww 

XXX 

WWW 

VVVV 

UUUUUUUU 

UUUUUUUU 

VVVV 

WWW 

X XX  Y 

YY  Y 

XXX 

wwww 

VVVVVVV 

VVVVVVV 

wwww 

XXX 

Y Y Y 

Y YY 

xxxx 

wwwwwwww 

xxxx 

Y YY 

YYYY  XXXXXXXXXXXXXXXXXXXXXXXXX  YYYY 


F igure 

33. 

Z configuration,  0^  = 

45°, 

RBT  = -0.393  (-0 

.393), 

AL  = 

1.0 

cm,  A0_  = 0.00428  rad. 

The 

curvature  along 

the 

x-coordinate  has  been  reduced  to  approximately  zero. 
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Figure  34.  Z configuration,  0^  = 45°,  RBT  = 0 (-0.361),  AL  = 1.0  cm, 
A©2  = -0.00580  rad.  The  curvature  along  the  y-coordinate  has  been 
reduced  to  approximately  zero. 
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Figure  35.  Z configuration,  0^  - 45°,  RBT  = -0.361  (-0.361), 
AL  = 1.0  cm,  A0£  = -0.00580  rad.  The  curvature  along  the 
y-coordinate  has  been  reduced  to  approximately  zero. 
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Figure  36. 


U configuration,  0^ 


45°,  RBT  = 0 (1.102),  AL  = A62 
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Figure  37. 


U configuration,  6^  = 45c 


RBT  = 1.102  (1.10), 


AL  = A02  = 0. 
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Figure  38.  Z configuration,  6^  = 70°,  RBT  = 0 (-1.997). 
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Figure  39.  Z configuration,  0 = 70°,  RBT  = -0.799  (-1.997), 
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Figure  40.  Z configuration,  0-  = 70°,  RBT  = 0 (-2.166),  A0  = +0.007. 
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Figure  41.  Z configuration,  0 = 70°,  RBT  = 0 (-1.837),  AO  = -0.007. 
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Figure  42. 


U configuration,  0^  = 70°,  RBT  = 0 (3.698), 
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Figure  43.  U configuration,  0^  = 70°,  RBT  = 1.48  (3.698), 


AL  = A62  = 0. 
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Figure  44.  Predicted  values  (smooth  curves)  and  calculated  values  of  the  root-mean-square  OPD 
in  microns  as  a function  of  AL,  for  Z configuration,  9 = 45°,  A9„  = 0. 
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e 46.  Predicted  values  (smooth  curves)  and  calculated  values  of  the  root-mean-square  OPD 
microns  as  a function  of  AL,  for  Z configuration,  0.  = 55°,  A0  = 0. 
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Figure  47.  Experimental  arrangement  for  obtaining  inter ferograms . 
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NO  THEORETICAL  OPD  AVAILABLE 
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Figure  51(b).  Z configuration,  0^  = 45°,  RST  = -0.296  (-0.370) 
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Z configuration 


Figure  53(b).  Z configuration,  0 


t*  • . 


Figure  55(a).  Misalignment  in  L,  &L 
9-1  = 45°,  Z configuration. 
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Z configuration,  0 


Figure  56(a).  Misalignment  in  AQ^  = 7 mrad, 
6^  = 45°,  Z configuration. 
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Figure  56(b).  0 = 45°,  RBT  = 0 (-0.392), 

A02  = 0.007  rad,  Z conf iguration . 
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Figure  57(a).  Misalignment  in  0„  + A0^  = -7  mrad , 
0 = 45°,  Z configuration. 
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Figure  57(b).  Z configuration,  9 = 45°, 

RBT  = 0 (-0.348),  A62  = -0.007  rad. 
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Figure  61.  Variation  of  RBT , =70°,  Z configuration. 


NO  THEORETICAL  OPD  AVAILABLE 


Figure  62(a).  Variation  of  RBT 


Figure  62(b).  RBT  = 0 (-1.997),  0 


Figure  64.  Variation  of  RBT,  6^  = 70°,  Z configuration. 
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NO  THEORETICAL  OPD  AVAILABLE 
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Z configuration 


Figure  66.  Misalignment  in  I 


NO  THEORETICAL  OPD  AVAILABLE 
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Figure  70.  Misalignment  in  0£>  A0^  “ 14  mraa, 


0 = 70  , Z configuration. 


NO  THEORETICAL  OPD  AVAILABLE 


U configuration 
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BT  = 0 (3.698). 


Figure  72(a).  Variation  of  RBT,  0^  = 70°,  U configuration. 


Figure  72(b).  U configuration,  0^  = 70°,  RBT  = 1.48  (3.698) 
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Figure  74.  Tilt  in  6^  in  a plane  orthogonal  to  the  plane  of  rotation 
of  and  02>  AB^  = -7  mrad,  B^  = 70°,  Z configuration. 


NO  THEORETICAL  OPD  AVAILABLE 
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Appendix  A.  EQUATIONS  FOR  THREE-DIMENSIONAL  RAY  TRACING 


The  coordinate  system  used  at  a particular  stage  of  ray-tracing 
calculations  will  ordinarily  be  chosen  such  that  the  central  ray  is 
along  the  positive  z axis.  Usually  the  origin  of  coordinates  is  at  the 
intersection  of  the  central  ray  with  the  spherical-mirror  surface.  The 
equations  can  be  related  to  Figure  A-l,  which  depicts  the  case  of  a 
convex  mirror.  Provision  is  made  for  the  mirror  to  be  tilted  about  the 
x-axis,  i.e.,  in  the  y-z  plane,  by  an  angle  T.  The  coordinate  system 
is  left-handed. 


The  basic  ray-optics  equation  for  reflection  is  simply 


e = e . 
r 1 


2 n cos  9 , 


(A-l) 


where  e^  and  e are  unit  vectors  along  the  incident  and  reflected  rays, 


respectively,  n is  a unit  normal  vector  (pointed  from  the  intersection 
point  toward  the  mirror  center) , and  6 is  the  angle  between  the  incident 


(or  reflected)  ray  and  the  mirror  normal  n. 


The  mirror  tilt  angle  T is  taken  as  positive  for  counterclockwise 
rotation,  i.e,.,  from  the  nositive  z direction  toward  the  positive 
y direction.  The  convention  is  adopted  that  the  mirror  curvature  R is 
negative  for  convex  mirrors.  The  coordinates  of  the  mirror  center  are 
then  0,  -R  sin  T,  -R  cos  T;  this  is  valid  for  either  sign  of  R.  The 
equation  of  the  mirror  surface  is 


2 2 2 2 
(x,y,z)  = const  = R =x  +(y+R  sin  T)  + (z  + R cos  T) 


(A-2) 


The  unit  normal  vector  can  be  determined  by  taking  the  gradient  of  this 
function,  i.e., 


S ‘ S « ’ 5 ’ l + sl"  T-  f + “s  T 


(A- 3) 


Now  z can  be  eliminated: 


= -(f  + > /-(f)2 T) 


(A-4) 


The  direction  cosines  of  the  incident  ray  are  denoted  by  a,  b,  and  c, 
i.e.,  giving  vector  equation: 


et  = (a,b,c) 


(A-5) 
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and  similarly,  (using  primes  to  denote  the  reflected  ray) , 


er  = (a' ,b' ,c')  . (A-6) 

The  cosine  of  the  angle  0 between  the  incident  ray  and  the  surface 
normal  is  given  by  the  dot  product  of  the  corresponding  unit  vectors, 
that  is: 

cos  0 = e.  . n . (A-7) 

1 

This  gives  the  explicit  result: 

cos  = 0 a (f)+  b(f)  + c / 1 -(f)2  - (*  + sin  t)2  . (A-8) 

This  equation  can  conveniently  be  used  to  determine  the  numerical  value 
of  9.  The  direction  cosines  of  the  reflected  ray  can  then  be  determined 
from  the  following  equations: 

a'  = a - 2 cos  0 • ^ , (A-9) 

(A-10) 


+ sin  T^ 2 • (A-ll) 

For  the  special  case  of  an  incident  ray  parallel  to  the  z-axis, 
i.e.,  for  a = b = 0,  the  following  is  noted: 


c'  = c - 2 cos2  0 . (A-12) 

The  necessary  relationship  between  the  direction  cosines 

(a')2  + (b')2  + (c')2  = 1 (A-13) 

is  reserved  here  for  a possible  test,  but  is  not  used  to  determine  the 
value  of  c1,  because  an  Equation  (11-A)  which  determines  the  sign  as 
well  as  the  magnitude  of  c1  is  needed. 

Equations  (9-A)  through  (11-A)  for  direction  cosines  of  reflected 
ray  require  knowledge  not  only  of  direction  cosines  of  incident  ray  but 
also  of  the  coordinates  of  the  point  of  intersection  of  the  incident  ray 
with  the  mirro^r.  In  general,  these  must  be  determined  from  a separate 
calculation.  For  the  special  case  of  the  incident  ray  parallel  to  the 


z-axis,  encountered  in  starting  some  ray-tracing  problems,  x and  y will 
be  given,  and  z can  be  determined  from  the  equation  obtained  by  solving 
Equation  (A-2)  for  z: 


z = -R  ccs  T - 


2 2 
x - (y  + R sin  T) 


(A-14) 


At  various  stages  in  the  calculations,  it  will  be  desirable  to  make 
rotational  coordinate  transformation.  In  particular,  such  a rotation 
of  coordinates  will  be  made  after  each  reflection  from  a spherical 
mirror.  As  indicated  in  Figure  A-2,  the  direction  of  the  central  ray  is 
changed  such  that  the  associated  angle  in  the  z-y  plane  is  increased  by 
it  + 2 T,  where  T is  the  mirror  tilt  angle,  taken  as  positive  for  coun- 
terclockwise rotations.  This  is  denoted  by  e- , i.e.:  • 


= TT  + 2 T 


(A-15) 


With  reference  to  Figure  A-3,  a rotational  transformation  from  an 
unprimed  to  a primed  coordinate  system  as  shown  has  the  transformation 
equations : 


z'  = z cos  6'  + y sin  ^ , (A-16) 

y'  = -z  sin  ■&  + y cos  & . (A-17) 


The  effect  of  the  coordinate  rotation  on  direction  cosines  can  be 
derived  by  reference  to  Figure  A-4.  The  polar  and  azimuthal  angles 
V and  W define  the  direction  of  the  ray  vector  g,  which  is  alternatively 
defined  by  the  direction  cosines  a,  b,  and  c.  Thus  the  following  can 


be  written: 

a = cos  V , (A-18) 

b = sin  V sin  W , (A-19) 

c = sin  V cos  W . (A-20) 

These  equations  can  be  readily  solved  for  V and  W: 

V = cos_1(a)  , (A-21) 

" ■ sin'1  (nb)  ■ c°rl  (nr*)  • (A-22) 


The  effect  of  a rotation  about  the  x-axis,  measured  as  positive  in  the 
sense  discussed  previously  is  to  reduce  W algebraically  by  the  magnitude 
of  the  rotation  angle,  while  the  polar  angle  V is  unaffected.  Thus,  the 
following  can  be  written: 


Figure  A-l.  Convex  mirror 


The  direction  cosines  can  then  be  written  in  the  rotated  frame  in  terms 
of  the  new  angles: 


a"  = sin(V)  = a 
b"  = sin(V)  sin(W') 
c"  = sin(V)  cos(W') 


(A-25) 

(A-26) 

(A-27) 


Here  double  primes  are  used  on  the  a,  b,  and  c simply  to  distinguish 
from  the  equations  relating  to  the  change  in  direction  cosines  as  a 
result  of  reflection  from  a spherical  mirror. 

For  computer  evaluation  of  W,  the  proper  quadrant  must  be  chosen. 

A suitable  method  is  evaluating  the  magnitude  of  W from  the  second 

part  of  Equation  (A-22)  involving  the  cos  ^ function,  but  with  the  sign 
of  W chosen  to  conform  to  the  first  part  of  Equation  (A-22) , involving 

the  sin  ^ function. 

After  finding  the  direction  cosines  and  coordinates  at  the  mirror- 
intersection  point,  referred  to  a coordinate  system  with  positive  z-axis 
along  central  ray,  the  case  that  the  next  mirror  is  also  a spherical 
mirror  is  considered.  First  the  point  of  intersection  of  the  ray  with 
this  mirror  is  found.  For  this  purpose,  the  parametric  equations  with 
t are  written  for  the  ray  as: 


x = x + a t 
y'  = y + b t 
z'  = z + c t 


(A-28) 

(A-29) 

(A-30) 


Equation  (A-2)  can  be  rewritten  defining  the  mirror  surface: 

(x’)2  + (y'  + R sin  T)2  + (z’  + R cos  T)2  = R2 
Substituting  Equations  (A-28)  through  (A-30)  gives: 


(A-31) 


(x  + at)2  + (y  + R sin  T + bt)2  + (z  + R cos  T + ct)2  = R2. 
This  can  be  expressed  as: 


(A-32) 


t + 2Bt  + C = 0 


(A-33) 


where 


B = ax  + b(y  + R sin  T)  + c(z  + R cos  T)  , 
C = x2  + (y  + R sin  T)2  + (z  + R cos  T)2  - R2 


(A- 34) 
(A-35) 


The  solution  to  Equation  (A-33)  is  simply: 


The  choice  of  sign  in  front  of  the  square  root  should  be  such  as  to 
obtain  a positive  solution.  It  should  be  the  smaller  in  magnitude  of 
the  two  in  the  event  that  there  are  two  positive  roots. 


Appendix  B.  EFFECTS  OF  INCREMENTS  TO  MIRROR  SPACING  L 
AND  TO  CONCAVE-MIRROR  TILT  ANGLE  fl2 


By  the  nature  of  the  concept,  there  will  be  no  net  focusing  or 
defocusing,  i.e.,  no  second-power  term  in  x or  y in  the  equation  for  a 
constant-phase  surface,  in  the  output  when  both  L and  0^  are  set  pre- 
cisely at  their  nominal  values.  This  fact  is  of  considerable  practical 
importance.  But  because  there  may  be  small  departures  of  AL  and  A02 

from  their  nominal  values,  it  is  appropriate  to  investigate  the  pre- 
dicted effects  (i.e.,  predicted  in  a fairly  simple  closed-form  analysis 
as  given  in  this  appendix)  of  such  departures  and  to  compare  these 
predictions  to  results  of  detailed  calculations. 

First,  increments  AL  are  considered.  The  effective  focal  lengths 
of  the  mirrors  in  the  x-z  plane  are  denoted,  by  F^  and  F,^.  Similar 

notation  will  be  employed  for  the  y-z  plane;  differences  in  the  two 
cases  arise  only  from  the  presence  of  cos  0^  or  sec  0^.  F denotes 

the  effective  focal  length  in  the  x-z  plane  of  the  output  beam.  The 
following  lens-type  equation  can  then  be  written: 

1 


FTx  Flx  - L - 4L  F2x 

F2x  + (F1»  - L - tL>  . 

<Flx  - L - SL)F2x 

But  L has  a value  such  that  the  output  is  collimated,  i.e., 

F2x  + Flx  " L = °* 

It  is  here  understood  that  F^^  is  negative.  Then,  to  a good  approxima- 
tion, the  following  can  be  written: 

1 AL 

F ~ 2 

T*  <F2x> 

Now  the  incremental  optical  path  length  for  the  case  of  a wave  focusing 
with  focal  length  F.^,  is  simply 

z = C2Q  x2 


Similarly,  for  focusing  in  the  y-z  plane,  the  following  is  obtained: 


The  corresponding  value  of  C^q  is  given  by 


The  specific  expression  for  is 


F2x  ‘ T ,ec  *2 


hence, 

C 


2AL  2 0 
20--R^2COS  9 2 


Similarly,  the  following  is  obtained: 

_ 2AL  2 

^02  * “ 2 sec  ®*>  * 

R2 

A numerical  example  with  R2  = 675  cm,  9^  = 55°,  0^  = 38.76°  has  been 

considered.  The  predictions  of  the  preceding  equations  for  this  example 
are 


I 


C2Q  = -0.267  x 10”5  AL  , 

CQ2  = -0.722  x 10-5  AL 

The  values  of  C2q  and  C^2  have  also  been  determined  by  simple  function- 
fitting to  results  of  numerical  OPD  calculations.  A comparison  of  these 
latter  calculated  results  with  the  predicted  results  derived  is  given 
in  Table  B-l.  The  agreement  is  satisfactory. 


» 


-r-  Jl 
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It  should  be  noted  that  the  signs  of  and  CQ2  are  the  same  for 

increments  in  mirror-spacing  L.  That  is  to  say,  the  net  focusing  or 
defocusing  (converging  or  diverging)  effect  is  qualitatively  the  same 
(though  quantitatively  different)  in  the  two  transverse  planes. 


Next,  the  effect  of  an  incremental  tilt  angle  A@2  is  considered. 
The  equation  for  focal  length  of  output  in  x-z  plane  is  now 

1 1 1 


Tx 


F.  - L 
lx 


F„  + AF„ 
2x  2x 


which  reduces  to  a high  degree  of  approximation  to 


Tx 


AF2x 

‘V2 


but. 


2x 


~2  sec  02 


and  hence 


AF2x  = — sec  02  tan  02  A62 


This  leads  to 


~rh~  = _ sin  e2  a02 

FTx  R2  1 


but  again, 
C20  = 


2F, 


Tx 


This  finally  gives 


C = — sin  6„  A0_ 
°20  R2  22 
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For  the  yz  plane,  the  following  is  obtained: 


F2y T C°S  92 

The  following  is  readily  obtained: 

R 

F2y  = T (_sin  02)  A62 

This  finally  gives 

c02  = " ^ Sln  92  S6c2  (62)  A02 

In  the  preceding  equations  A0^  is  assumed  to  be  in  radians. 

A numerical  example  with  = -290  cm,  R2  = 675  cm, 

0^  = 0.95993  rad  (55°),  and  02  = 0.67647  rad  (38.76°)  has  been  consid- 
ered. The  predicted  values  of  the  quadratic  parameters  C2q  and  C02  are 

C20  = 0.927  X lO-3  A@2 

CQ2  = -1.525  x 10-3  A62 

A comparison  of  values  predicted  in  this  way  with  values  obtained  by 
function-fitting  to  numerically  calculated  OPDs  is  given  in  Table  B-2. 
The  agreement  is  satisfactory. 

C2q  and  CQ2  have  opposite  signs  for  the  case  of  increments  in  tilt 

angle  02>  This  is  readily  understood  by  noting  that  the  effective  focal 

lengths  in  the  two  transverse  planes  are  changed  in  opposite  directions 
by  such  an  incremental  tilt  angle.  This  behavior  is  in  contrast  to  the 
behavior  with  regard  to  increments  in  mirror  spacing. 
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TABLE  B-l.  COMPARISON  OF  PREDICTED  AND  CALCULATED  VALUES  OF  BEAM 
CURVATURE  PARAMETERS  C2Q  AND  CQ2  RESULTING  FROM  INCREMENTS  AL  IN 

MIRROR  SEPARATION.  (Rx  = -290  cm,  R2  = 675  cm,  6^^  = 55°, 

©2  = 38.76°,  L = 180.01005  cm)* 


*Predicted  values  are  the  same  for  Z and  U configurations;  calcu- 
lated values  are  for  Z configurations. 


TABLE  B-2.  COMPARISON  OF  PREDICTED  AND  CALCULATED  VALUES  OF  BEAM 
CURVATURE  PARAMETERS  C2Q  AND  CQ2  RESULTING  FROM  INCREMENTS  A62 

IN  CONCAVE-MIRROR  TILT  ANGLE.  (RL  = -290  cm,  R£  = 675  cm, 

6 * 0.95993  rad  (55°),  02  = 0.67647  rad  (38.76°), 

L = 180.01005  cm)* 
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TABLE  B-l . COMPARISON  OF  PREDICTED  AND  CALCULATED  VALUES  OF  BEAM 
CURVATURE  PARAMETERS  C AND  CQ2  RESULTING  FROM  INCREMENTS  AL  IN 

MIRROR  SEPARATION.  (R1  = -290  cm,  R2  = 675  cm,  0^^  = 55°, 

02  = 38.76°,  L = 180.01005  cm)* 


L (cm) 

C20 

C 

20 

Predicted 

Calculated 

Predicted 

Calculated 

+1.0 

-0.267E-5 

-0.268E-5 

-0.722E-5 

-0.728E-5 

-1.0 

+0.267E-5 

+0.266E-5 

+0.722E-5 

+0.716E-5 

+2.0 

-0.534E-5 

-0.536E-5 

-1.444E-5 

-1.457E-5 

-2.0 

+0.534E-5 

+0.531E-5 

+1.444E-5 

+1.430E-5 

*Predicted  values  are  the  same  for  Z and  U configurations;  calcu- 
lated values  are  for  Z configurations. 


TABLE  B-2 . COMPARISON  OF  PREDICTED  AND  CALCULATED  VALUES  OF  BEAM 
CURVATURE  PARAMETERS  C2Q  AND  CQ2  RESULTING  FROM  INCREMENTS  A02 

IN  CONCAVE-MIRROR  TILT  ANGLE.  (Rj  = -290  cm,  R£  = 675  cm, 

= 0.95993  rad  (55°),  02  = 0.67647  rad  (38.76°), 

L = 180.01005  cm)* 


02 

C20 

C02 

Predicted 

Calculated 

Predicted 

Calculated 

+0.0035 

+0.324E-5 

+0.325E-5 

-0.534E-5 

-0.540E-5 

-0.0035 

-0.324E-5 

-0.324E-5 

+0.534E-5 

+0.528E-5  \ 

+0.0070 

+0.648E-5 

+0.652E-5 

-1.068E-5 

-0.082 E-5 

-0.0070 

-0.648E-5 

-0.646E-5 

+1.068E-5 

+1.054E-5 

*Predicted  values  are  the  same  for  Z and  U configurations;  calcu- 
lated values  are  for  Z configurations.  Note  that  0.0035  rad 
0.2°. 
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